Implementing bvp4c in MATLAB